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We examine three different ways a heavy quarkonium can dissociate at temperatures below the quark-gluon plasma phase transition 
temperature T c : spontaneous dissociation, dissociation by thermalization, and dissociation by collision with hadrons. We evaluate 
the cross sections for the dissociation of J/'ifj and T in collision with n as a function of temperature, using the quark-interchange model 
of Barnes and Swanson. We also evaluate the dissociation temperatures for various quarkonia, and the fraction of quarkonium lying 
above the dissociation threshold as a function of temperature, using an interquark potential inferred from lattice gauge calculations. 


I. Introduction 

The suppression of heavy quarkonium production has 
been the subject of intense interest as it .was proposed 
as a signature for the quark-gluon plasmatH. The recent 
experimental observation of an anomalous J/ip suppre&j 
sion in Pb+Pb collisions by the NASO Collaboration □ 
has been considered by many authors 0. However, there 
is considerable uncertainty on the origin of the anoma¬ 
lous suppression due to the lack of reliable experimental 
and theoretical information on J/ip and \J dissociation 
below the quark-gluon plasma phase transition tempera¬ 
ture T c . 

We shall study three quarkonium dissociation mech¬ 
anisms at T < T c . A heavy quarkonium can dissociate 
spontaneously when it becomes unbound at a tempera¬ 
ture above its dissociation temperature. A quarkonium 
system in thermal equilibrium with the medium can dis¬ 
sociate by thermalization when a fraction of the quarko¬ 
nium lies at energies above the dissociation threshold. A 
heavy quarkonium can also dissociate by collision with 
light hadrons. We shall discuss these three mechanisms 
in turn. Some of the details of the discussions can be 
found in Refs. [4] and [5]. 

II. Dissociation Cross Section at T=0 

For heavy quarkonium at T = 0, the only dissociation 
mechanism for J/ip, x> and T is through their colli¬ 
sion with hadrons. To date, three general approaches 
have been proposed for the calculation jpf the dissoci¬ 
ation cross sections: perturbative QCD 0, effective La- 
grangiansQ, and a quark-intojjphange model using explicit 
quark-model wavefunctions o@. The cross sections esti¬ 
mated using these approaches vary over many orders of 
magnitude. 


We favor the use of the quarkpinterchange model pro¬ 
posed by Barnes and Swanson □ , E2l, as it has previously 
been successfully applied to many analogous light-quark 
systems Ed where it is in good agreement with experi¬ 
ment. The generalization of the model to finite temper¬ 
atures can also be carried out when the wave function 
of the heavy|jquarkonium is determined as a function of 
temperatures. A relativistic generalization of the quark- 
based model has also been presentedEix 

To obtain the dissociation cross sections, we first use 
the known hadron masses to determine interquark forces 
and bound state wave functions. We write the interaction 
in the form 


V(r) 


A(i) A (j) f a s _ 3 
2 ' 2 { r 4 r 


87ra s 

3mimj 



+ Vcc 


(1) 


where for an antiquark, the generator A/2 is replaced by 
— \ T /2. From the meson mass spectrum, we obtain the 
following set of parameters: 

12tt 

a " ~ (33 — 2n/) ln(A + Q 2 /B 2 ) ’ ’ 

A = 10, B = 0.31 GeV, Q 2 = (bound state mass) 2 , 
b = 0.18 GeV 2 , cr = 0.897 GeV, V con = 0.620 GeV, 
m u = m d = 0.334 GeV, m s = 0.575 GeV, 
m c = 1.776 GeV, and m b = 5.102 GeV. 


The dissociation of a heavy quarkonium in collision 
with a meson can be described as a quark-interchange 
process. The Born-order scattering amplitude can be 
evaluated as overlap integrals of hadron bound state 
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wavefunctions using the “Feynman rules” given in Ap¬ 
pendix C of Ref. [9]. 



collision with 7r, K, and p at T =0 as a func¬ 
tion of the kinetic energy Eke- 

We can calculate the scattering amplitude in the 
“prior” formalism in which the interaction takes place 
before the interchange of the quarks. There is the corre¬ 
sponding “post” formalism in which the interaction takes 
place after the interchange of the quarks. In nonrelativis- 
tic bound state scattering, one finds that these results 
are actually equal, provided that the same Hamiltonian 
is used to generate lie asymptotic bound states and to 
drive the scattering td. We will adopt an intermediate 
approach and assume relativistic kinematics and phase 
space but use nonrelativistic scattering amplitudes; in 
consequence, we find different ‘post’ and ‘prior’ cross sec¬ 
tions in general. Here we will take the mean value of the 
‘post’ and ‘prior’ results as our estimated cross section; 
separate ‘post’ and ‘prior’ cross sections will be shown in 
Figs. 1-3 as an indication of our systematic uncertainty. 

When the interaction Eq. (jlj) is used to calculate the 
low-energy reaction, 1 = 2 tttt scattering phase shifts, ex¬ 
perimental data agree well with the theoretical model. In 
addition the quark-interchange model has been applied 
successfully to KN and NN scattering, as well as to a 
wide rang£.of related no-annihilation K*N , AIV, and AN 
reactionsEJ. 

Dissociation cross sections at T = 0 calculated with 
the Barnes and Swanson model are shown in Figs. 1-3. 


Fig. 1 is for J/tjj, Fig. 2 for % c2 , and Fig. 3 for T. 



Fig. 2 Dissociation cross sections of \c 2 in col¬ 
lision with 7r, I \, and p at T =0 as a function 
of the kinetic energy Eke- 



Fig. 3 Dissociation cross sections of T in col¬ 
lision with 7r, K , and p at T =0 as a function 
of the kinetic energy Eke- 
One finds from Figs. 1-3 that the cross sections are 
sensitive functions of energy. The maxima of the cross 
sections for tt + K + J/ip, and 7r + \c 2 are about 
1-2 mb, and the cross sections for exothermic reactions 
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p + J/’ip, K + Xc2, and p + Xc2 decreases from about 6 
mb to about 0-2 mb as Eke decrease from 0.05 GeV to 
0.2 GeV. The maxima of the cross sections of 7r, K, and 
ponT are about 0.1 mb. 


III. Heavy Quarkonium at T <T C 

Recently, Digal, Petreczky, and Satz 0 reported theo¬ 
retical results on the dissociation temperatures of heavy 
quarkonia. The basic input is the temperature depen¬ 
dence of the Q-Q interaction as inferred from lattice 
gauge calculations of Karsch et alx3. 

Due to the action of dynamical quark pairs, a proper 
description of the heavy quarkonium state sh ould be 
based on a screening potential even at T = 0 as 
the heavy quarkonium becomes a pair of open charm or 
open bottom mesons when r becomes very large. 

To study the behavior of a heavy quarkonium at 
finite temperatures, we calculate the energy e of the 
heavy quarkonium singlc^particle state ( QQ)jls from 
the Schrodinger equation El 

f h 2 

j-V- V + V (r, T) + (toq + mg 

- Mq 9 - M q q)6(R - r)^ JLS {r) = eip JLS (r). (3) 


The two-body system consists of a quark pair Q and Q 
at a distance R and becomes a pair of lowest-mass 
mesons Qq and qQ at r > R. The energy e is measured 
relative to the pair of lowest-mass mesons at r —> oo. 
The quarkonium is bound if e is negative. It is unbound 
if e is positive and dissociates spontaneously, subject to 
selection rules. The quantity p \2 is the reduced mass. 
We use itiq, uiq , Mq q , and rn q Q at T = 0 and include all 
the temperature dependence of the effective interaction 
in V (r, T). For numerical calculations, R has been set to 
°.S fm. 

Following Karsch et al. liJ, we represent the Q-Q- ifc 
teraction by a Yukawa plus an exponential potential! 1 url 


V(r,T ) 


W^ (r)r _ b(T) (T)r 
3 r p(T) 


(4) 


where b(T ) is the effective string-tension coefficient, p(T) 
is the effective screening parameter, and the potential is 
calibrated to vanish as r approaches infinity-. The results 
of the lattice calculations of Karsch et al .lid can be de¬ 
scribed by p(T) = po = 0.28 GeV and 

b(T) = b 0 [l-(T/T c ) 2 ]9(T c -T) : (5) 

where bo = 0.35 GeV 2 . The value of po has been fixed 
to be the same as at T = 0, and the value of bo is close 
to the b value obtained earlier at T = Olid. 

In our calculation, we use a running a s as in Eq. (j^). 
We include, in addition, the spin-spin, spin-orbit, and 


tensor interactions, with a running spin-spin interaction 
width as in Ref. [18]. The eigenvalues of the Hamilto¬ 
nian can be obtained by matrix diagonalization using a 
set of nonorthogonal Gaussian basis states with different 
widths, as described in Refs. [10,4,5]. 

IV. Spontaneous Dissociation 

Figure 4 shows the charmonium single-particle states as 
a function of temperature 0 . In considering the sponta¬ 
neous dissociation below T c , it is necessary to find the 
selection rules for the dissociation of a heavy quarko¬ 
nium state with initial quantum numbers J, Lj, and S t 
into two mesons with a total spin S and a relative orbital 
angular momentum L. With interactions of the type in 
Eq. (® and the spin-spin interaction, the total S is con¬ 
served. Parity conservation requires A L = \L — Lj| = 1 
in this dissociation. Hence, J/ip and xj)' will dissociate 
into a pair of open charm mesons with L = 1, while x 
states will dissociate into a pair of open charm mesons 
with L = 0 or 2. These conservation laws give rise to 
selection rules for the spontaneous dissociation of heavy 
quarkoniumEI. The difference in the energy of the lowest 
allowed final state and the energy of the pair of lowest- 
mass meson state gives the additional increment of the 
threshold energy. 



Fig. 4. Charmonium single-particle states as 
a function of temperature. The threshold en¬ 
ergies are indicated as horizontal lines. The 
solid circles indicate the locations of the dis¬ 
sociation temperatures. 

The dissociation temperatures of relevant quarko¬ 
nium states can be determined by plotting the state en¬ 
ergies of the quarkonia and the threshold energies after 
considering the selection rules. The points of intercept, 
as indicated by solid circles in Fig. 4, give the positions of 
the dissociation temperatures. Figure 5 shows the bot- 
tomium single-particle states as a function of T/T c . 
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We list the dissociation temperatures of charmonia 
and bottomia in Tables I and II. With the additional 
threshold energies due to the angular momentum selec¬ 
tion rules, the dissociation temperatures are shifted to 
higher energies. The increase in the dissociation temper¬ 
ature is large for charmonia and small for bottomia. 

We confirm the general features of the results of Di- 
gal et al. but there are also differences as the dissocia¬ 
tion temperatures depend on the potential and the selec¬ 
tion rules. With the temperature-dependent potential of 
Eqs. ©-©. the dissociation temperatures of all heavy 
quarkonia, except \bo , Xbi> Xb 2 , and T, are below T c , 



Fig. 5. Bottomium single-particle energies as 
a function of temperature. The solid circles 
indicate the locations of the dissociation tem¬ 
peratures. 


Table I. Charmonium dissociation tempera¬ 
tures Td in units of T c 


Charmonium 

V 

Xc2 

Xcl 

J/iP 

T d /T c 

0.50 

0.91 

0.90 

0.99 

T d /T c 

(Selection Rules 
not invoked) 

0.37 

0.66 

0.72 

0.94 

T d /T c 

(Digal et al.) 

0.1-0.2 

0.74 

0.74 

1.10 


Table II. Bottomium dissociation temperatures 
Td in units of T c 


Bottomium 

'Y'// 

Xb2 

Xm 

Y'/ 

Xb2 

X61 

T 

TdlT c 

0.57 

0.82 

0.82 

0.96 

>1.0 

>1.0 

>1.0 

TdlT c 

(Selection Rules 
not invoked) 

0.54 

0.75 

0.78 

0.95 

1.00 

>1.00 

>1.0 

Td/Tc 

(Digal et al.) 

0.75 

0.83 

0.83 

1.10 

1.13 

1.13 

2.31 


V. Dissociation by Thermalization 

If the heavy quarkonium is in thermal equilibrium with 
the medium, the occupation probabilities of the heavy 
quarkonium state e ? ; will be distributed according to the 
Bose-Einstein distribution, 

exp{(ej — aO/T} — 1 ’ ^ 

where /i is the chemical potential. Such a thermal 
equilibrium arises from inelastic reactions of the type 
h+(QQ)jLS —> h'+(QQ)j>L'S'- When the heavy quarko¬ 
nium reaches thermal equilibrium with the medium, 
there is a finite fraction of the quarkonium system to 
be found in excited states. The fraction / of the system 
lying above their thresholds for spontaneous dissociation 
will dissociate into open charm or open bottom mesons. 
We call such a process dissociation by thermalization. 



Fig. 6. The fractions of charmonium and bot¬ 
tomium lying above the dissociation threshold 
as a function of T/T c . 


We evaluate the fraction / as a function of temper¬ 
ature for charmonium and bottomium, and the results 
are shown in Fig. 6. A state label along the curves de¬ 
notes the onset of the occurrence when that state emerges 
above its dissociation threshold. As one observes, the 
fraction / increases with temperature. As the temper¬ 
ature approaches T c , the fraction of charmonium lying 
above dissociation thresholds is larger than the corre¬ 
sponding fraction of bottomium. 

VI. Dissociation by Collision with Pions 

We study in this section the dissociation oi J/ip and T 
by collision with pions in a medium at temperature T. 
As the temperature of the medium increases, the quarko¬ 
nium single-particle energy changes and the dissociation 
threshold energy decreases. As a consequence, the disso¬ 
ciation cross section will change. 
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We can calculate the dissociation cross sections of 
J/ip and T in collision with 7 r as a function o£the temper¬ 
ature using the Barnes and Swanson model 0. The calcu¬ 
lation requires the energies and wave functions of the ini¬ 
tial and final states, as well as the interquark interaction 
which leads to the dissociation. For the interquark inter¬ 
action, we generalize the temperature-dependent Yukawa 
and exponential interaction in Eqs. ([|) and (j|) by replac¬ 
ing the color factor —4/3 there with the color operator 
(\(i)/2) • (X(j)/2) where for an antiquark, the generator 
A/2 is replaced by —A T /2. 

Using the single-particle energies and wave func¬ 
tions obtained with this temperature-dependent interac¬ 
tion, the scattering amplitude can be evaluated and the 
dissociation cross section can be calculated, as in Ref. 
[4].The sum of dissociation cross sections for tt + J/ip —> 
DD *, D*D , D*D* are shown in Fig. 7a for different tem¬ 
peratures T/T c , as a function of the kinetic energy Eke- 
In Fig. 7b we show similar total dissociation cross sections 
for 7 r + T —> BB*, B*B, B*B* for various temperatures. 



< GeV ) 

Fig. 7. Total dissociation cross section oi J/ip 
(Fig. a) and T (Fig. b) in collision with n for 
various temperatures as a function of the ki¬ 
netic energy Eke- 

We observe in Fig. 7 that the maximum values of the 
dissociation cross sections increase and the positions of 
the maxima shift to lower kinetic energies as the tem¬ 
perature increases. Such an increase arises from the de¬ 
crease of the threshold energies as the temperature in¬ 
creases. Over a large range of temperatures below the 
phase transition temperature, dissociation cross sections 
of J/ip and T in collisions with n are large. 

In a hadron gas, pions collide with a heavy quarko- 
nium at different energies. We can get an idea of the 
energy-averaged magnitude of the dissociation cross sec¬ 


tion by treating the pions as a Bose-Einstein gas at tem¬ 
perature T. In Fig. 9, we show the quantity (av) which 
is the product of the dissociation cross section and the 
relative velocity averaged over the energies of the pions. 
The quantity (av) is about 2 mb at T/T c = 0.70 and 
rises to about 3 mb at T/T c = 0.95 where the value of T c 
has been taken to be 0.175 GeVEj. 



Fig. 8. The average of the dissociation cross 
section of J/ip and T in collision with i r mul¬ 
tiplied by v as a function of the temperature. 


We can estimate the survival probability of a heavy 
quarkonium in a hot pion gas in the presence of this 
type of collisional dissociation. If we represent the sur¬ 
vival probability S by exp{—/}, the exponential factor I 
is given by I = J^ freeze (av)(r)p(r)dT, where p(r) is the 
density of ir at the proper time r, and tq and Tf ree ze are 
the initial proper time and the freeze-out proper time re¬ 
spectively. The quantity (av) in Fig. 8 can be represented 
approximately by 

(av) = (av) c (T/T c ), (7) 


where (av) c ~ 3 mb for J/ip or T. Assuming a Bjorken 
expansion, we obtain 


I = 


(&v)c 


dN 

dyA 



f Pfreeze^O A \ 

V dN/dy J 


( 8 ) 


where Tq is the initial temperature and pf reeze is the 
freeze-out pion density. 

The degree of absorption by collision with pions de¬ 
pends on the initial absorption time To, the freeze-out 
pion density pf ree ze) and the initial temperature Tq. If 
To ~ T c , pfreeze = 0.5/ fm 3 , and To = 3 fm/c+R/ 7 , then 
for the most central Pb-Pb collision at 158A GeV, the 
heavy quarkonium survival probability is S ~ 0.5 and 
for the most central Au-Au collision RHIC at y' snn = 
200 GeV, the heavy quarkonium survival probability is 
S ~ 0.1. These estimates show that the absorption of 
both J/ip and T by the hot pion gas is substantial. 
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VII. Discussions and Conclusions 


We study the dissociation of a heavy quarkonium at T < 
T c . At T = 0, dissociation is possible only by collision 
with hadrons. We calculate the dissociation cross section 
for the dissociation of J/ij), Xi and Y i n collision with n, 
K, and p at T = 0. The cross sections for J/■i/j and T 
dissociation in collision with 7 r are relatively small. 

The temperature of the medium alters the interac¬ 
tion between a heavy quark Q and antiquark Q in the 
medium. We calculate the quarkonium single-particle 
states using a potential inferred from lattice gauge cal¬ 
culations of Karsch et all3 and obtain the dissociation 
temperatures. We confirm the general features of the re¬ 
sults of Digal et ablifl but there are differences. We find 
that the selection rules change the dissociation temper¬ 
atures substantially for charmonia but only slightly for 
bottomia, and that the dissociation temperatures of all 
heavy quarkonia, except Xi> 0 ) X6i> Xb 2 , and T, are below 
T c . 

A quarkonium in a medium can collide with par¬ 
ticles in the medium to reach thermal equilibrium. A 
heavy quarkonium in thermal equilibrium can dissociate 
by thermalization as there is a finite probability for the 
system to be in an excited state lying above its dissocia¬ 
tion threshold. We find that the fraction of the quarko¬ 
nium lying above the threshold increases with increasing 
temperatures. 

Dissociation of a heavy quarkonium can occur in col¬ 
lision with hadrons. As the temperature increases, the 
threshold energies for collisional dissociation decrease. 
As a consequence, the dissociation cross sections increase. 
We have estimated the absorption of J/ if> and T in col¬ 
lision with pions in central Pb-Pb collisions at SPS and 
RHIC energies and found the absorption to be substan¬ 
tial. Further microscopic investigations of the dissoci¬ 
ation of heavy quarkonium in collision with hadrons in 
high-energy heavy-ion collisions will be of great interest. 
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